Evolution of the bursting- layer wave during a Type 1 X-ray burst 
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ABSTRACT 

In a popular scenario due to Heyl, quasi periodic oscillations (QPOs) which are seen during type 1 X-ray 
bursts are produced by giant travelling waves in neutron-star oceans. Piro and Bildsten have proposed that 
during the burst cooling the wave in the bursting layer may convert into a deep crustal interface wave, which 
would cut off the visible QPOs. This cut-off would help explain the magnitude of the QPO frequency drift, 
which is otherwise overpredicted by a factor of several in Heyl's scenario. In this paper, we study the coupling 
between the bursting layer and the deep ocean. The coupling turns out to be weak and only a small fraction 
of the surface- wave energy gets transferred to that of the crustal-interface wave during the burst. Thus the 
crustal-interface wave plays no dynamical role during the burst, and no early QPO cut-off should occur. 

Subject headings: Type 1 X-ray burst, Quasi periodic oscillations, shallow-water wave, crustal-interface wave 



1. Introduction 

Unstable nuclear explosions of accreted hydrogen and 
helium on the surface of a neutron star (NS) produce ob- 
servable (type-1) X-ray bursts in some low-mass x-ray bi- 
naries. The bursts typically have rise times of order 1 s 
and decay times of order 10 — 100 s. During some of the 
bursts quasi-periodic oscillations (QPOs) are seen with fre- 
quencies ranging from 40 to 700 Hz. In the decaying tail 
of the burst the QPO frequency typically shows a positive 
drift of 1 — 5 Hz to an asymptotic value, which is char- 
acteristic to a bursting object 0. This stability has led 
many to conclude that this asymptotic frequency is equal 
to the spin of the NS; in this picture the QPO is pro- 
duced by a rotating asymmetric temperature pattern on 
the NS surface. However, after more than 10 years since 
the discovery of the burst QPOs, it is still unclear what 
leads to the surface asymmetry at late times on these low 



1 Notable exceptions are burst QPOs from millisecond X-ray pulsars. 
The QPO frequency of these objects stays constant during the bursts 
and is equal to the NS spin. 



magnetic-field NSs. The frequency drifts are also challeng- 
ing to explain, since the NS spin should remain constant 
during the burst. Strohmayer et al. 1997b argued that the 
conservation of angular momentum could be the cause of 
the drifts. In their scenario the spin of the bursting layer 
changes due to radial expansion and contraction during 
the burst. However detailed calculations show that this 
effect under-predicts the frequency drift in the burst tail 
by a factor of 2-3 (Cumming et al. 2002). It also does not 
explain the surface asymmetry at late times when the hot 
spot has already spread over the entire surface. 

In an important paper Heyl (2004) proposed the non- 
radial surface modes (waves) as the origin of the QPOs. In 
his scenario the burst excites waves in the bursting layer, 
and it is these waves that modulate the X-ray flux period- 
ically. The QPO's frequency drift can then be accounted 
for by the change of the wave frequency during the cooling 
of the bursting layer. However, the basic surface wave fre- 
quencies change by order 10 Hz during the burst cooling, 
significantly larger than the observed amount. Two pro- 
posals have been put forward to overcome this difficulty. 
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Firstly, Heyl (2004) has proposed that the modes localised 
in the photosphere will have smaller frequencies and drifts. 
However, due to the small density and vertical scaleheight, 
the photosphcric mode frequencies are extremely sensitive 
to the presence of hydromagnetic stresses in the neutron- 
star ocean. For example, even for a small vertical field of 
B = 10 5 G, the Alfven crossing frequency 

v a ~ B ~ 30(B/10 5 G)(10cm/fe)[10 4 g/cm/p] 1/2 Hz, 
zhy/Anp 

(1) 

is much greater than the frequency of the original mode. 
Here p and h are the characteristic density and scaleheight 
of the photosphere. Thus even a small magnetic field is 
expected to strongly modify the photospheric modes and 
substantially increase their frequencies. 

Secondly, Piro and Bildsten (2005b) (hereafter P&B05b) 
have suggested that the global surface wave is responsible 
for the QPOs, but that the frequency drift is cut short by 
resonant conversion of the bursting-layer wave (BLW) to 
the crustal-interface wave (CIW). The latter is an oscilla- 
tory mode whose energy is concentrated near the interface 
between the crystallised (crust) and liquid (deep ocean) 
burst ashes (Piro & Bildsten 2005a [P&B05a]). The CIW 
has a low frequency of ~ 4Hz, and its main restoring force 
is due to the shear modulus of the crystallised ashes. The 
frequency of the CIW is unperturbed during the burst, 
which implies a stable asymptotic frequency of the QPOs. 
One of the interesting features of this scenario is that the 
QPO asymptotic frequency differs from the NS rotational 
frequency by that of the CIW. This would have important 
implications for the narrow-band LIGO searches for the 
NS-generated gravitational waves. 

In this paper we investigate in detail the P&B05b sce- 
nario. First, we evaluate the coupling between the bursting 
layer and the deep ocean in a shallow-water model. We find 
a good qualitative agreement with modal P&B05b calcula- 
tions for a more realistic ocean model. Second, we analyse 
in detail the dynamics of the mode conversion. P&B05b 
have conjectured, without proof, that the cooling period 
is long enough for the resonant wave conversion to occur. 
By contrast, we show that the cooling period is in fact too 
short for the resonant wave conversion to occur, and thus 
that P&B05b burst-QPO scenario is not viable. 

2. The bursting- layer wave and the crustal- 
interface wave 

In Hcyl's (2004) scenario a type 1 X-ray burst excites 
a backwards-propagating buoyant R-mode in the bursting 
layer. During the cooling period the pressure scale height 



of the bursting layer decreases. This makes the frequency 
of the wave decrease, which will be seen as a positive drift 
by an inertial observer. The deep ocean remains unper- 
turbed during the burst so the frequency of the CIW is un- 
changed. At a certain point during burst cooling the wave 
frequencies will approach each other. However the two 
waves represent coupled vibrational modes of the NS. Be- 
cause of the coupling, the vibrational mode frequencies do 
not cross but undergo an avoided crossing (P&B05b). The 
mechanism for this phenomenon will be explained later on; 
for an illustration of an avoided crossing between modes see 
figured] At the avoided crossing the two modes exchange 
their characteristics: the mode that is the BLW becomes 
the CIW and vice verse. If the mode frequencies evolve 
sufficiently slowly (adiabatically) with time, the system, 
having started in one mode, remains continuously in the 
same mode. This implies that upon a sufficiently slow pas- 
sage through the avoided crossing, the BLW converts into 
CIW. P&B05b have argued that this is what happens dur- 
ing the burst cooling, and they have argued that the con- 
version cuts short the frequency drift of the BLW. While 
we agree with P&B05b about the location and strength 
of the avoided crossing, we show below that the evolution 
is strongly non-adiabatic through the crossing, and hence 
that the BLW^CIW conversion does not occur. 

The BLW rides on top of the burst ashes. Taking a thin- 
ocean limit (see Pedlosky [1979] )B the dispersion relation 
is 

u*~gH b k 2 ^, (2) 
P 

where g is the gravitational acceleration, Hj, is the pres- 
sure scale height at the interface between the bursting layer 
and the ashes, k is the effective wave number (to be de- 
fined later for the case of rapid rotation) , p is the density of 
the ashes at the interface, and Ap is the magnitude of the 
density jump across the interface. Here the Ap/p is impor- 
tant, since the density difference between the bottom of the 
bursting layer and top of the ocean becomes smaller during 
the burst cooling, which contributes to slowing down the 
wave. The radial component of the Coriolis force is much 
smaller than the buoyancy force (Bildsten, Ushomirsky & 
Cutler 1996, here after BUC96), which allows one to sepa- 
rate the angular and radial parts of the non-radial modes. 
Thus the mode solutions are as follows: 

V (r, 9, cb, t) = g(r)k(fi, 4>) exp(-iwt), (3) 

where r is the radius, 9 is the latitude, 4> is the longitude, 
t is the time, T)(r, 9, <f>, t) is the vibrational mode, g(r) is 

2 In this limit the radius of the NS is much bigger than the thickness 
of the bursting layer 
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the radial part of the mode, fc(0, (ft) is the angular part 
of the mode and u> is the angular frequency of the mode. 
Because the NS surface is azimuthally symmetric the az- 
imuthal part is periodic: 



r){r, 9, <j>, t) = g{r)f{9) exp i(m(f> - ut), 



(4) 



where f(9) is the latitudinal part of the mode and m the 
azimuthal eigenvalue. The latitudinal part is given by the 
Laplace's tidal equation: 



J q.m 



[/ g ,m(0)] = -A(<z,m)/,, m (0) 



(5) 



with L the Laplace tidal operator, q — 2il/u> the ratio be- 
tween the NS rotation and the angular frequency of the 
mode, 9 the latitude, / a latitudinal eigenfunction and A 
the corresponding eigenvalue. The quantity X/R 2 is the 
effective wavenumber k for the waves (Longuet-Higgins 
1968). The mathematical expression for the Laplace tidal 
operator is as follows, 



L„ 



d_ 



1-C 2 d 
1 - q 2 ( 2 d( 



(i-C 2 )(i-g 2 C 2 ) 

qm(l + q 2 C 2 ) 

2/-2\2 ' 



(6) 



(i-<z 2 C 2 ) 



where ( — cos 9. Longuet-Higgins (1968) found the angu- 
lar eigenfunctions, called Hough functions, in terms of the 
spherical-harmonics series. Later BUC96 gave a simpler 
method to find the eigenfunctions numerically. To approx- 
imately match the QPO with the rotational frequency, the 
angular part should be that of the m = 1 buoyant R-mode 
(m = 1, I = 2 R-mode for slow rotation; see Heyl 2004). 
This mode has A ~ 0.11 in the fast rotating limit (oj 17; 
BUC96). P&B05b estimated the initial frequency of the 
BLW to be around 10.8 Hz in the frame rotating with the 
NS, with a negative drift of about 9 Hz during the cooling 
period H. 

The energy of the CIW is concentrated near the crustal 
interface with the deep ocean. The dispersion relation is 



ul,. 



gH c k 2 



(J) 



where [i is the shear modulus of the crust and P and H c 
are the pressure and pressure scale height, respectively, at 



3 To find these values they considered the bursting layer as an 
ideal gas, whose temperature is set by a vertical radiative flux of 
10 25 ergs.cm~ 2 .s _1 . This layer extends till a depth where the den- 
sity equals 10 6 g/cm 3 and the pressure scale height is about 2 m. 
During the burst cooling the temperature of the layer changes from 
order 10 9 K to order 10 8 K. 



the crustal interface (P&B05a). The factor /i/P is present, 
since the crustal elasticity provides the restoring force for 
the CIW. Crystallisation of the electron degenerate ocean 
material is determined by the dimensionless parameter 



L = 



(Zef 
aksT 



(8) 



where fcf, is Boltzmann's constant, T is the temperature, 
a = (3/47rrii) 2 , rii is the ion number density, e is the elec- 
tron charge and Z is the average ion charge. Following 
Farouki & Hamaguchi 1993, the crystallisation sets in when 
r ~ 173, at a density of about 10 9 g/cm 3 . The pres- 
sure scale height at the crustal interface is of order 20 m. 
P&B05b have obtained a frequency of 4.3 Hz for the CIW^. 

3. The shallow-water model 

P&B05b have computed the mode frequencies of the 
coupled BLW and CIW and have demonstrated the exis- 
tence of the avoided crossing. In their computation, they 
have assumed a certain 1-dimensional model for the ver- 
tical density profile of the bursting layer. However, this 
profile is poorly known in detail. The deflagration during 
the burst is a multi-dimensional phenomenon, with most 
of the fuel ignited from the top by a strongly inclined defla- 
gration front (Spitkovsky, Levin, & Ushomirsky 2002), and 
thus any 1-dimensional model is probably inaccurate. In 
this section we work out a simple shallow-water model for 
the coupling between layers, and find the strength of the 
avoided crossing, and hence the coupling strength between 
the modes. Our results are in good qualitative agreement 
with P&B05b; thus the existence and the strength of the 
avoided crossing is a generic feature of the bursting ocean, 
and is insensitive to the uncertain details of the model. The 
rest of this section is rather technical, and can be safely 
skipped by the reader not interested in the mathematical 
details. 

Our basic idea is to represent the bursting layer and 
deep ocean by two layers of incompressible inviscid fluid 
(figure [l}. The bottom beneath the two layers is flexible, 
and represents the crust. The top of the bursting layer is 
at the radius R + h t , the top of the deep ocean is at R + h m 
and the crustal interface is at R+hf, (hb = in equilibrium 
state). For this model we make use of the shallow-ocean 
approximation. In this situation the thickness of the fluid 
layers, which is of order H, is small compared to the radius 
R of the NS, i.e. we assume S = H/R <C 1. The dynam- 
ical equations of motion for the thin water layers are the 



4 The wave frequency was calculated for a NS accreting He-rich mate- 
rial at an accretion rate M > W~ 9 M<7). 
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1 Bursting layer 




Fig. 1. — This image shows a cut out of the shallow-water 
model. Two layers on a large spinning star representing 
the bursting layer and the deep ocean of NS surface. The 
bursting layer has height h t and the deep ocean has height 
h m above the crust. 



continuity equation for an incompressible fluid 



and the momentum equation for an inviscid fluid 

dv - ^ -VP 

— + 20 x v = V$, 

at p 



(9) 



(10) 



where v is the velocity and <£> is the gravitational potential. 
There is no exchange of fluid between the layers. Following 
Pedlosky (1979) and expanding Eq. (|10p to zeroth order in 
S we find the radius-independent equations for the trans- 
verse velocity: 



dvgi 
dt 

dVB2 

dt 



2Q,v^\ cos 9 = 
2f2u02 cos 9 = 



-gd s h t 

R 1 
-gd e h m gpi 



dv. 



4,1 



dt 
dv^ 
dt 



2flvgi COSt 

2f2vg2 cost 



R 

-gd^ht 

Rsm9 ' 
Rsin9 



dg(h t - h m ), 



gpid^jht - hm) 
Rp2 sin 9 



(11a) 

(lib) 
(11c) 

(lid) 



Here the labels 1 and 2 stand for the top and the bottom 
layers, respectively. To solve for the radial velocity we 
substitute Eqs. fTTajl . !|llb|) . ([TTc|) and (jlldj) in the conti- 
nuity equation @ and the boundary conditions mentioned 
above. We keep the terms up to the first order in 6, since 



the radial velocity is of this order relative to the horizontal 
one. We find 



Bi r , i , dh m 
v rl = -[h m -y} + — 

B2 r , i , dhb 



(12a) 
(12b) 



where y is the height above the bottom of the shallow- water 
model, so that r = R + y. From the continuity equation, 
we find 



dh t dh m _ Bi 

~dT~^r-li [hm ~ htU 



dh r 



dt 



dhb B2 r , 

'W = T [hb ~ hm]> 



(13a) 
(13b) 



where B is the angular part of the velocity divergence 



B = 



sin 9 



+ dgVg + COt 9vg 



(14) 



Now we work out the modes using a linear perturbation 
theory. We look for solution in the following form: 

Q(r,6,(f>,t) = Q (r) + Q'(r,9)expi[m<j)-LJt}, (15) 

where Q is a dynamical variable, Qq is its equilibrium 
value, and Q' is the Eulerian perturbation amplitude. Us- 
ing Eq. (fT5|) in the Eqs. (|llap and (|llc|) and keeping only 
linear terms we find 



i2ujQ cos( 



ui 2 ^\ + i2ujfl cos 9^gi — 



R ' 
-igmh' t 

RsmO ' 



(16a) 
(16b) 



where £ is the amplitude of the Lagrangian displacement 
vector. Using Eqs. (|16ap and (|16bp we can solve for £^1 and 
£01 in terms of h' t . Perturbing equation (fT4| and expressing 
£01 and £01 in terms of h' t , we get 



(17) 



Analogous to the derivation for layer 1 we find for layer 2 



B' = 



'IE. 
Rw J 



h' m + ^(h' t -h' m ) 

92 



From Eqs. (fT3"a|) and (|13b|l we find 

iu [h' m - h' t 



B' 

-p- [h m o — hto] 



iv [h' b - h' m 



R 



(18) 

(19a) 
(19b) 
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These differential equations can be separated in a radial 
part and an angular part. The angular part is Laplace's 
tidal equation (see eq. [5]). Rewriting Eqs. (|19a,|l and (|19bj) 
using Eqs. (fT7)l . (fTS|) and ([5]) we find the frequency of the 
modes. There are two equations 



g\(q, m) [h t0 - h m0 ] 



1 



R 2 



('-ft) 



g\(q, m) [h m0 - h b0 ] 






t 1 - 















and (20a) 



(20b) 



In the shallow-water model the depth within a layer is 
equal to the pressure scale height. Thus, as the term with 
(hto — h m o) indicates, Eq. (|20ap describes a wave concen- 



trated in the bursting layer, i.e. the BLW. The vertical 
displacement of the top of the ashes affects the frequency 
of this wave, as is represented in Eq. (|20ap by the term 
h' m /h' t . Eq. (|20b[) describes a wave concentrated in the 



deep ocean, i.e. the CIW. The flexing of the crust is rep- 
resented by the term h' b /h' t in Eq. (|20b|) . 

3.1. Finding the avoided crossing 

During a type 1 X-ray burst the bursting layer expands 
and a density discontinuity is created at the bursting layer 
and ocean interface. This density discontinuity produces a 
buoyancy force at this interface, which acts as the restoring 
force for the BLW. The dispersion relation for the BLW is 
given by Eq. $Z§, The equation contains a factor Ap/p, 
which mimics the strength of the buoyancy at the top of 
the ashes. The standard dispersion relation for a shallow- 
water wave is 

to 2 = g Hk 2 , (21) 

which applies to waves in the shallow-water model. This 
equation does not contain the factor Ap/p of Eq. @. To 
let the model take this term into account we define the 
effective pressure scale height of the BLW to be: 



Hi 



A P 



(22) 



In the shallow- water model we set the pressure scale height 
of the bursting layer to H b r , rather than H^. During the 
burst cooling on the NS the pressure scale height of the 
bursting layer and the density discontinuity at the bursting 
layer and ocean interface decrease. To imitate both these 
effects in the model we simply decrease the value of the 
pressure scale height of the top layer. 

We use the following representative parameters for our 
shallow-water model. 



1. The angular frequency of the NS is fi ~ 1000 rad.s . 

2. On the surface of the NS the centripetal force is just 
one percent of the strength of the gravitational force. 
We approximate the radial acceleration as the grav- 
itational acceleration; g ~ 1.858 x 10 14 cm.s~ 2 . 

3. The mass of the bursting layer is 2.55 x 10 21 g. 

4. The mass of the deep ocean is 2.52 x 10 25 g. 

5. The pressure scale height of the bursting layer is set 
to 2.28m to give the BLW a frequency of 10.8 Hz. 

6. The pressure scale height of the deep ocean is set to 
40.59 m to give the CIW a frequency of 4.3 Hz. 

7. The angular eigenvalue X(q, m) is set to that of the 
m=l buoyant R-mode. 

For the calculations of the masses of the layers, we used 
an equation of state with P ~ p 4 / 3 . This applies to the 
burning layer, which consists of an ideal gas and where 
the flux is approximately constant within the layer, and to 
the deep ocean, which has a pressure dominated by rela- 
tivistic degenerate electrons. The equation of hydrostatic 
equilibrium is 

dP , . 

* = ~ 9P (23) 
With this equation and the equation of state of above the 
density can be rewritten as a function of the radius 



p~a(R~r)\ 



(24) 



where a is some constant. We solved for a using the NS 
surface parameters given by P&B05b, and by integrating 
over the density we found the masses of the layers. To 
find the latitudinal eigenvalue A for the m = 1 buoyant 
R-mode, we solved 'Laplace Tidal equation' using the nu- 
merical method described in BUC96. 

A pressure perturbation in the deep ocean causes the 
neutron star crust to flex. To incorporate this flexing into 
our model we followed the analysis of P&B05a. We as- 
sumed that the ocean is plane-parallel instead of spherical 
and we denote the horizontal and vertical coordinates by 
x and z, respectively. The bottom of the model represents 
the crustal interface of the NS surface. The zz component 
of the Lagrangian perturbation of the stress tensor Aa zz , 
must be continuous at the crustal interface to avoid infinite 
vertical acceleration. On the crustal side of the interface, 
we have 

3ri/4 + ^ p2j (25) 



Aa° 



1 



hto V 
Send / 
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where Ti is the adiabatic exponent and s enc i is the depth 
from the surface where the CIW ends (e(s en d) = 0). On 
the fluid side of the crustal interface, 



Act? 



bove 



= -pig{K - Ki) - P2g(K n - h' b ). 



(26) 



Equalising equations ([2"ij]) and pi?)) , setting r ~ 4/3 and 
taking s en d >> hto (see P&B05a) we find the following 
relation for the bottom flexing 



hi 



Pi 



hL 



_ I 1 m m 



(27) 



We express h' b /h' t in terms of h' m /h' t using Eq. (|27p and 
solve Eqs. (|20ap and (|20b|l for the frequencies of the modes 
and h' m /h' t . 
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Fig. 2. — In this figure the frequencies of the 2 modes of the 
model have been plotted versus the frequency of the uncou- 
pled BLW mode, lui/(2it). The modes undergo an avoided 
crossing. At the avoided crossing mode 1 changes from the 
BLW to the CIW and mode 2 vice verse. 

In figure [2] the frequencies of the two modes have been 
plotted versus the , uncoupled BLW mode, Qx/(2tt), for a 
NS spinning at 1000 rad/s. We determine the presence of 
an avoided crossing, and measure the minimal frequency 
difference (i.e., the gap) between the modes to be 0.042Hz. 
This is similar, in order of magnitude, to that in P&B05b, 
who have measured 0.058Hz for a more realistic model of 
the NS ocean. 

In the next chapter we will use the size in the gap of the 
avoided crossing to find the coupling between the bursting 
layer and deep ocean, and to model the dynamics of the 
system as it passes through the avoided crossing. 



4. Dynamics of transition through the avoided 
crossing 

The avoided crossing found in the previous section re- 
sults from the coupling between the top and the bottom 
layers. Without the coupling, each layer would have a vi- 
brational mode frequency set by the individual dispersion 
relation for shallow- water waves. When coupling is intro- 
duced the layers push and pull each other. However, away 
from the avoided crossing, the coupling is weak and the 
modes are identified with the individual layers. P&B05b 
argued that as H c g is slowly reduced, at the avoided cross- 
ing the two modes would adiabatically exchange their char- 
acteristics. The mode concentrated in the bursting layer 
would transfer its energy to the deep ocean and vice verse. 
P&B05b have used the following criterion for the adia- 
baticity: r S> 1/u), where r is the characteristic timescale 
of the frequency change. This criterion is fine away from 
the crossing; however, as is well-known from quantum me- 
chanics, it breaks down near the avoided crossing. The 
theory of passing through the avoided crossing was worked 
out independently by Landau (1932) and Zener (1932). In 
this section we first model explicitly the mode crossing nu- 
merically, and then show that our results are consistent 
with the Landau-Zener theory. Our conclusions are oppo- 
site to those reached by P&B05b: the burst cooling is not 
slow enough for the modes to adiabatically exchange their 
characteristics. 

Let X\ and x 2 be the generalised vibrational amplitudes 
of the two modes localised in the bursting layer and deep 
ocean respectively. These modes are coupled via a general 
interaction Hamiltonian i?i nt = — ax\x 2 , and their joint 
equations of motion are 



m\x\ + k\X\ = ax2, 
m 2 x 2 + k 2 x 2 — axi, 



(28a) 
(28b) 



where mi is the effective mass of the mode in the bursting 
layer, m 2 is the effective mass of the mode in the deep 
ocean, fej/m, is the vibrational frequency of the layers set 
by the dispersion relation for shallow-water waves. It is 
convenient to rewrite these equations as follows: 



xi + u x xi = ax 2 , 
x 2 + lj 2 x 2 — axi, 



where Xi — Xiy/rrii, a = aj ^jm\m 2 and u>i = ki/nii 
two proper frequencies of the system are given by 



(29a) 
(29b) 

The 



± 



(30) 
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Note that 

\{J+-J-) = j& + {Ql-<4) 2 . (31) 
Therefore, the effective coupling a can be expressed as 

(32) 



Energy transfer between the oscillators 



a 

7^2 



1 (col - u 2 _ 



Thus can be read off directly from the plot in Fig. [5J 

(33) 



Auj 



where Aw is the minimal separation between the two 
branches. Using Eq. (|33|) and the data of Fig.[2]we find the 
coupling between the layers for different rotational spin of 
the NS which is shown in table [1] Compared to the angu- 



Spin (rad.s- 1 ) (0) 


coupling (rad.s '*) (a) 


1000 


7.29 


4000 


7.41 



Table 1: This table shows the effective coupling coefficient 
between the bursting layer and deep ocean. 

lar frequencies of the BLW and CIW the coupling is rather 
small, i.e. ^ « 1. This is because of the huge effective 
mass difference between the two layers: the deep ocean is 
~ 10 times heavier than the bursting layer. 

We now integrate Eqs. (|29ap and (|29bp numerically, us- 
ing the values of a in table [T] We set the initial frequencies 
of oscillator 1 and 2 equal to the frequencies of the burst- 
ing layer wave (10.8 Hz) and the crustal interface wave 
(4.3 Hz), respectively. To simulate the burst cooling we let 
the frequency of oscillator 1 decrease exponentially by 9 Hz 
during a time period set by the typically observed cooling- 
period of 10 s. More precisely, we evolve the frequency of 
oscillator 1 as a function of time as follows: first, we keep 
it constant for 10 seconds, then we evolve it for 10 seconds 
according to the following expression 



2ir 



14.24 Hz exp 



-{t- 10s) 
10s 



3.44 Hz, (34) 



and then we keep it constant for another 10 seconds. In 
figure[3]we plot the energy of the two oscillators normalised 
to the initial energy of oscillator 1, over the whole of the 
30 s time interval. Oscillator 1 starts with a certain am- 
plitude, while oscillator 2 begins with zero displacement 
and velocity. Because initially the coupling between the 
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Fig. 3. — This figure shows the energy transfer between 
two coupled harmonic oscillators. During the first and last 
10 s the frequencies of the oscillators are constant, but in 
the range 10 — 20 s the frequency of oscillator 1 goes down 
crossing the frequency of oscillator 2. 



oscillators is small both oscillators are close to vibrational 
eigenmodes of the system, with oscillator 1 representing 
the excited BLW at the start of the burst cooling period. 
During the first 10 s virtually no energy is transferred be- 
tween the oscillators. In the second interval of 10 s the 
energy of oscillator 1 decreases and the energy of oscilla- 
tor 2 increases, especially at the location of the avoided 
crossing. However the energy of oscillator 2 remains much 
smaller than the energy of oscillator 1. In the last 10 s 
again there is virtually no energy transfer. We conclude 
that during the burst cooling the system does not evolve 
adiabatically and that the BLW does not change into a 
CIW. 

We have run the simulation for several different values 
of yito, d and mi/m?, and in each case found that the 
energy transfer between the oscillators is very small. The 
most important parameter is the cooling time; with cooling 
times > 100 s we get a significant decrease for the BLW 
energy. However, this is far longer than most observed 
cooling times. Thus mode conversion cannot be a robust 
feature of the Type 1 X-ray bursts. 

We now compare our results with the Landau-Zener the- 
ory (Landau 1932, Zener 1932; for a quick derivation, see 
Wittig 2005). In the Landau-Zener problem the Hamilto- 
nian is changing so that the difference between the energy 
levels of a quantum-mechanical system is decreasing lin- 
early with time and then crossing zero. This gives the sys- 
tem, which is originally in the first energy state, a chance 
to make a transition to the second energy state. The ex- 
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pression for the transition probability is 



dimming, A. et al. 2002, ApJ, 564, 343 

Farouki, R., & Hamaguchi, S. 1993, Phys. Rev. E, 47, 4330 
Heyl, J.S. 2004, ApJ, 600, 939 
Landau, L.D. 1932, Phys. Z., 2, 46 
Lee, U. 2004, ApJ, 600, 914 



( 



d(u)i - uo 2 ) 



) 



1 



pi^ 2 = exp 




dt 



(35) 



Using the formula we find that p\^2 is close to unity for our Longuet-Higgins, M. S. 1968, Phil. Trans. R. Soc. London, 262, 511 

parameters, and thus that resonant conversion is not ex- Pedlosky, J. 1979, Geophysical Fluid Dynamics (New york:Springer) 

pected to occur. The conversion would occur if pi^2 were Piro, A. L. , & Bildsten, L. 2005a, ApJ, 619, 1054 (P&B05a) 

to remain close to zero; for this the cooling time greater Piro, A. L. , & Bildsten, L. 2005b, ApJ, 629, 438 (P&B05b) 

than ~ 100s is required, in agreement with our numerical Spitkovsky, A., Levin, Y., & Ushomirsky, G. 2002, ApJ, 566, 101 

calculations. Strohmayer, T. E. et al. 1991, ApJ, 375, 679 



Wittig, C. 2005, J. Phys. Chem. B, 109, 8428 

In this paper we have shown that the bursting-layer Zenerj c 1932j Proc . R . So c. London A, 137, 696 
wave does not turn into the crustal-interface wave dur- 
ing Type 1 X-ray burst cooling. This argues against the 
P&B05b picture in which the asymptotic frequency of the 
burst QPO differs from the NS spin by the crustal-interface 
wave frequency (4 Hz) . Our conclusions are consistent with 
a very recent discovery that the accreting neutron star Aql 
X-l is a transient millisecond x-ray pulsar (Casella et al., 
2007). In that object, the asyptotic maximal frequency of 
the burst QPOs is very close (~ 0.5Hz) to the spin fre- 
quency of the neutron star, without any shift due to the 
CIW. We draw two astrophysical conclusions from our cal- 
culations. Firstly, the simple surface waves do not work as 
the source of the burst QPO's. To find the source one 
should appeal to the waves in the photosphere (Heyl 2004) 
or some other non-trivial modes. Secondly, the asymp- 
totic spin of the QPOs is likely to be close to the NS spin, 
without modification due to the CIW motion. Advanced 
LIGO should take this into account for the future searches 
of gravitational waves from accreting and bursting NSs. 
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